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KINETIC  EQUATIONS  FOR  PLASMAS1 
by 

'i' 

Toyoki  Koga 

Polytechnic  Institute  of  Brooklyn,  Graduate  Center 
Farmingdale,  New  York 

SUMMARY 

Because  of  the  time  scale  assumption  necessary 
for  its  derivation,  the  Bogoliubov-Born-Green-Kirkwood- 
Yvon  (BBGKY)  hierarchy  of  equations  is  not  applicable  to 
systems  consisting  of  charged  particles,  except  for  special 
subsystems.  A  new  class  of  equations  governing  the 
evolutions  of  charged  particles  is  derived  from  the  Liouville 
equation  and  is  coarsed-grained  with  respect  to  time  and 
similar  particles.  In  the  zeroth  approximation,  there  are 
two  basic  types  of  inter -particle  interactions:  One  is  of  the 
Vlasov  type  and  the  other  of  tie  Boltzmann  type  characterizing 
interactions  among  nearest  neighboring  particles.  In  higher 
order  approximations,  mutual  perturbations  among  those 
basic  interactions  result  in  secondary  effects;  for  example, 
two  nearest-neighboring  particles  exert  a  force  of  microscopic 
order  to  another  particle.  Depending  on  the  ratio  between  the 

^This  research  was  supported  under  Contract  Nonr  339(38)  for 
PRO  JECT  DEFENDER,  and  the  Advanced  Research  Projects 
Agency  under  Order  No.  529  through  the  Office  of  Naval  Research. 
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number  of  electrons  and  the  number  of  ions  in  a  real  system,  the 
simulating  model  varies.  The  main  purpose  of  the  paper 
is  to  present  schemes  of  rational  treatment,  rather  than 
to  provide  numerical  results  in  detail  for  a  particular 
system. 


I.  INTRODUCTION 


The  Boltzmann  collision  integral,  based  on  the 
binary  collision  assumption  and  the  Stosszahlansatz  (  the 
assumption  of  collision  number),  diverges  in  the  case  of 
charged  particles.  As  a  more  radical  approach,  many 
authors ^  have  attempted  to  investigate  ionized  gases  by 
treating  the  Bogoliubov-Born-Green-Kirkwood-Yvon  (BBGKY) 

hierarchy  of  equations  derived  from  the  Liouville  equation. 

2 

As  has  been  argued  often  by  the  present  author  ,  however,  the 
BBGKY  hierarchy  is  not  feasible  for  consideration  as  the  basis 
of  treating  ionized  gases  in  general,  because  of  the  time  scale 
assumption  necessary  for  the  derivation  of  the  hierarchy;  the 
time  scale  relation  necessary  for  the  derivation  differs,  except 
for  special  cases,  from  the  time  *cale  relation  conceivable 
in  a  real  system  of  charged  particles.  The  purpose  of  the 
present  article  is  to  propose  schemes  of  treating  the  Liouville 
equation  feasible  for  ionized  gases. 

In  Section  II,  we  shall  integrate  partly  the  Liouville 
equation  to  derive  equations  of  evolution  of  distributions  of 
sub-systems.  After  investigating  possible  time  scales  of 
dynamical  processes  in  real  systems  in  Section  III,  we  shall  coarse- 
grain  these  equations  for  sub  systems  with  respect  to  time  and 
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similar  particles  so  that  they  will  yield  kinetic  equations  for 
plasmas  in  the  following  sections.  Each  of  the  resultant 
equations  his  its  limited  regime  of  feasibility;  the  feasibility 
depends  on  the  nature  and  the  condition  of  an  ionized  gas  under 
consideration.  Providing  numerical  results  in  detail  for  a 
particular  system  is  left  for  future  publication. 
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II.  EQUATIONS  OF  EVOLUTION  OF  SUB -SYSTEMS 


We  consider  a  system  consisting  of  N  particles  which 
are  represented  by  material  points  with  masses  m^m^i..,  mN* 
The  distribution  of  the  system  in  the  6N-dimensional  phase  space 

(N)  (n> 

is  denoted  by  D'  \  Function  D  satisfies  the  Liouville  equation 


where  p^  is  the  momentum  of  particle  i,  r^  the  position  vector 
of  the  same  particle,  and 

--3  *  -  3 

O  io  V  i  j 

i3  the  total  force  exerted  on  particle  the  external  force, 

[jp.  ■  the  force  exerted  by  particle  j.  Since  we  are  investigating 

(N) 

the  evolution  of  a  single  system,  D  ’  is  known  to  be  given  by 


(N)  N  * 

DlIN'  =  -  6  (X  -  X  (t)  ) 

i  =  1 

ajc 

where  X.  denotes  the  six-dimensional  vector  p.,  r.,  and  X.  (t) 
i  mil' 

represents  the  trajectory  of  particle  i  obtained  by  solving  the 
equations  of  motion  of  the  N  particles^.  The  equations  of  motion 
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constitute  the  characteristic  equations  of  Eq.  (2.  !).  Of  course, 

(N) 

D  given  abo\  5  satisfies  Eq.  (2.  1).  It  is  noted  that  X  is 

j}c 

an  ndependent  variable  while  is  a  function  of  time*  We 
mav  define  F^'  s  by 

F(1)(X.;t)  =  V  f  D(N)  tt  dX. 

1  J  it  i  J 

F(2,(X.X.;t)  =  V°  \  D(N)  tt  dX 
1  J  J  Mi.j  k 

etc. 

where  V  is  invariant  and  denotes  the  volume  of  the  space  within 
which  the  system  is  known  to  exist  throughout  the  period  of  time  of 
our  investigation.  Since  the  velocity  of  a  particle  is  finite.V  is 
also  finite.  It  is  convenient  to  define 


F(1)(i) 


F(1)(X.;t) 


F(a)(ij)  =  F{i)(X.  X.;t),  if  i  <  j 
=  F(i)(X.  X. ;t),  if  i  >  j 
=  0  ,  if  i  =  j 

etc. 
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(2.3) 


(2.3)' 


On  integration  of  each  term  of  Eq.  (2.  1)  with  respect  to  the 
phase-space  coordinates  of  all  the  particles,  except  for 
particle  i,  we  have 


f  i 


jL 

d  pi 


)  FC1) 


CD 


♦ 


z 


dX.,  ...dX,  =  0 


(2.4) 


Here  means  the  summation  of  similar  factors,  each 

of  which  rs  formulated  with  respect  to  one  of  the  M-particle 
subsystems  oade  by  ^rouping  the  particles,  except  for  particle  i, 
of  the  system.  Grouping  of  the  particles  to  v-particle  subsystems 
is  possible  only  when  (N-  1)/.  is  an  integer,  and  the  number  of  the 
subsystems  is  (N  -  l)/v  If  (N-1)/Vis  not  an  integer,  for  example 
101.  36.  we  have  to  admit  that  one  of  the  102  subsystems  consists 
of  less  than  panicles.  In  general  (N-  1)  »  v  ,  and  hence  the 
number  of  the  subsystems  is  large  so  that  we  may  ignore  the 
peculiarity  of  one  subsystem.  Therefore  we  always  treat  (N-l)/. 
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as  if  an  integer.  For  i,  there  are  N  different  choices.  When  i 
is  given,  the  number  of  different  ways  of  dividing  the  remaining  N  -  1 
particles  into  v  -particle  subsystems  is 


(N  -  i)i 


(N-l)/  , 

(v!)  L(N-l)v]i 


(2.5) 


In  the  microscopic  sense,  the  implication  of  Eq.  (2.  4)  does  not 
change  by  changing  the  way  of  grouping. 

It  is  easily  shown  that 
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f<>41) 


(ij...k)  dXJ...dXk 


1  i. 

rJV^i 


F(2)(ij)  dXj 


and  hence  we  may  write  for  Eq.  (2.  4) 


( 


6> 


c) 

*io’  **1 


♦ 


£ 

3pi 


F(2)(ij)  dX^ 


0 


(2.6) 


where  p<  tide  j  represents  all  the  particles  of  the  system  except 
for  particle  i.  As  is  well-known,  Eq.  6)  is  the  very  equation 
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which  yields,  after  being  coarse-grained,  the  first  equation  of 
the  BBGKY  hierarchy.  There  is  no  difference  between  Eq.  (2.  4) 
and  Eq.  (2.  6),  as  long  as  we  consider  them  in  the  precise  and 
microscopic  sense.  After  being  coarse-grained,  however,  the 
two  resultant  equations  are  different.  The  difference  is  due 
to  the  fact  that  the  equations  of  motion  of  particles  in  mutual 
.iteraction  are  non-linear,  s»nd  interactions  are  none  Multiple  tnan  binary* 
For  the  equation  of  evolution  of  F^1  wc  obtain, 
by  integrating  properly  Eq.  (2.  1) 


♦  •  •  •  • 


f(A1) 


(i J# • #k) 


♦ 


dX^...dXB  *  0 
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(2.6) 


In  the  above,  in  general,  it  is  not  trivial  to  choose  a 
proper  value  for  v'  •  As  has  been  stated  with  respect  to 
coarse -graining  of  Eq.  (2.4),  coarse -graining  of  Eq.  (2.6) 
with  a  different  value  for  leads  to  an  equation  of  different 
physical  implications.  We  assume,  however,  that  v  is 
sufficiently  large  that  the  effect  of  additional  v'  particles 
is  trivial;  this  assumption  is  of  a  nature  similar  to  that  of 
the  assumption  of  binary  collision  where  v  =  1  and  v'  =  0. 

Further  we  may  define  f  dX  by 

N  X  +dX 

f(X)dX  =  —  £  f  F  (1)(X.)dX.. 

V  11 

v  i=  1  X 

(2.7) 

f  dX  gives  the  number  of  the  particles  which  exist  in  the  states 
between  X  and  X  +  dX. 
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HI.  GENERAL  SCHEMES  OF  DERIVATION  OF  KINETIC  EQUATIONS 


First  of  all,  we  should  realize  that  the  limitations  of  our 
mathematical  ability  prevents  us  from  treating  precisely 
more  than  tv/o  body  problem  ..  Therefore  our  effort  must 
be  directed  toward  converting  N-body  problems  to  a  combination 
of  many  two-body  problems  in  an  approximation. 

Suppose  that  particle  1  is  the  test  particle  of  which 
the  time  evolution  of  our  main  concern.  One  of  the  field 
particles  is  named  by  2.  What  is  the  effect  of  particle  2  on 
the  test  particle?  If  particle  2  happens  to  come  close  to  the 
test  particle  within  a  certain  distance,  they  will  interact  with 
each  other  as  if  they  are  isolated  from  the  other  field  particles: 
their  mutual  interaction  is  much  stronger  than  their  interactions 
with  the  other  field  part'.cles.  But  such  a  strong  and  distinctive 
interaction  continues  to  exist  only  for  a  short  period  of  time. 
Particle  2  will  soon  be  remote  from  the  test  particle,  and  particle 
2  will  be  in  a  strong  interaction  with  its  new  nearest  neighbor, 
particle  3.  Of  course,  particle  2  does  yet  continue  to  exert  a 
force  on  the  test  particle.  But  the  behavior  of  particle  2 
is  almost  independent  of  the  behavior  of  particle  1  and  hence 
the  force  exerted  by  particle  2  on  particle  1  as  a  function  of  time. 


1  1 


changes  according  to  the  nature  of  the  interaction  between 
particle  2  and  particle  3;  the  behavior  of  particle  2  appears 
as  a  matter  of  fluctuation  from  the  view  point  of  particle  1. 

In  spite  of  such  fluctuations  appearing  in  the  behavior  of 
particle  2,  the  state  of  particle  2  is  steadily  localized  in 
a  larger  scale,  because  of  the  finiteness  of  its  velocity. 

Therefore,  the  average  of  the  fluctuating  force  due  to  particle 
2  over  a  proper  time  period  does  not  vanish;  particle  2  exerts 
a  weak  but  almost  stationary  force  on  particle  1.  Such  weak 
forces  due  to  many  field  particles  of  similarly  localized  states 
culminate  to  the  Vlasov  type  force  exerted  on  particle  1.  See  Fig.3»l» 
According  to  the  above  investigation,  there  are  at 
least  three  different  time  scales  of  the  force  exerted  on  particle  1 
by  particle  2:  1)  the  time  scale  of  the  close  interaction  be'veen 

the  test  particle  and  particle  2;  2)  the  time  scale  of  the  close  interaction 
between  a  field  particle  and  particle  2;  3)  the  time  scale  of  the 

average  fligh  of  particles  2.  If  all  the  panicles  are  of  the  same 
single  species,  time  scale  1  and  time  scale  2  are  of  the  same  order. 

On  the  other  hand,  if  the  system  under  consideration  is  composed 
of  particles  of  more  than  two  species,  the  time  scale  of  the  interaction 
between  the  test  particle  and  particle  2  may  be  different  from  the 
time  scale  of  the  interaction  between  particle  2  and  a  third  particle; 
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the  latte?  may  also  vary  depending  on  the  species  of  the  third 
particle.  Finally  the  time  scale  of  the  third  category  is 
almost  macroscopic. 

A  conclusion  derived  from  the  above  consideration  is 
that,  in  the  simplest  case  where  the  system  is  constituted  of 
particles  of  a  single  species,  it  is  necessary  to  take  for  v  in 
Eq.  (2.  4)  at  least  2  so  that  the  close  interaction  between  two 
neighboring  field  particles  is  taken  into  consideration  with  the 
same  time  scale  as  ti.at  of  the  close  interaction  between  the  test 
particle  and  a  field  particle.  This  caution  is  essential  for  coarse 
graining  Eq.  {2.  4)  with  respect  to  time. 

Since  it  has  been  realized  that  consideration  of  the  time 
scales  of  close  interactions  between  two  nearest  neighboring 
particles  is  essential,  we  investigate  the  orders  of  such  time 
scales  with  respect  to  typical  particles  in  the  following: 

By  assuming  that  the  gas  under  consideration  is  fully 

ionized  the  constituent  particles  are  ions  and  electrons.  It  is  also 

assumed  that  those  ions  are  of  a  single  species.  The  charge  of 

an  electron  is  denoted  by  -e  the  mass  by  m  ,  and  the  total  number 

e 

of  the  electrons  in  the  system  by  Ng.  The  charge  of  an  ion  is 
denoted  by  g  ,  the  mass  by  mj  and  the  total  number  of  ions  by 
N  j  .  In  general  we  may  assume  that 

'Ng  e  +  Nj  r  =  0  (3.  1) 

*j/  me  2000  (3.2) 
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We  also  assume  that 

e 

Due  to  relation  (3.  i)  the  velocity  of  an  electron 
much  larger  than  the  velocity  of  an  ion, not  only 
in  thermal  equilibrium,  but  also  when  the  state 
In  thermal  equilibrium,  it  may  be  easily  shown 


(3.3) 


is  usually 
when  the  gas  is 
is  non-uniform: 
that 


(the  average  magnitude  of  the  velocity  of  electron)  ^ 

(the  average  magnitude  of  the  velocity  of  an  ion) 

(3.4) 


If  the  gas  is  not  in  thermal  equilibrium,  the  situation  is  not  as 

simple.  However,  if  the  acceleration  of  an  electron  and  the 

acceleration  of  an  ion  occur  due  to  the  same  external  electro  — 

magnetic  field  relation  (3.  4)  may  be  a  reasonable  expectation, 

because  the  acceleration  of  an  ion  is  (r/e)(m  /m  )  times  larger 

®  I 

than  the  acceleration  of  an  electron. 

The  time  scale  of  the  interaction  between  two  nearest 
neighbroing  electrons  is  of  the  order  of 

-1/3 

7  =  n  /  «-c  n 

ee  e  e 

(3.5) 

where  n  is  the  number  density  of  electrons  and  <c  ">  is  the 
e  e 

av  age  speed  of  an  electron.  As  compared  with  the  average 
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speed  of  an  electron,  the  average  speed  of  an  ion  is  much 
smaller;  and  an  ion  appears  as  if  it  is  at  rest,  while  an  electron 
appears  to  move  swiftly.  The  time  scale  of  the  interaction 
between  an  electron  and  its  nearest  neighboring  ion  is  of  the 
order  of 


Tel  =  n  I  '  ~  e 
where  <  ce>  is  the  average  velocity  of  an  electron  and  n^ 

is  the  number  density  of  ions.  Since  n*and  n  are  assumed  to 

*  e 

be  of  the  same  order,  we  have 


1 

T* 


/  <  c  „> 


(3.6) 


T  e  I  —  Tee 


Cn  the  other  hand,  the  time  scale  of  the  interaction  between 
two  nearest-neighboring  ions  is  of  the  order  of 


1 

Ty 

T II  '  B  j  •  «  C  j  > 


(3,7) 


and  is  much  longer  fhan  t  ej  and  /  or  t  due  to  <<bt  %  <  c  > 


ee 


In  summation,  we  have 


b>>TsI-  T  ee 


(3.8) 
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IV.  MODEL  I 


4.  1.  General  Consideration 

The  first  model  is  a  system  constituted  of  similar 

particles  and  a  uniform  background  of  electric  charge  which 

compensates  the  total  charge  of  the  particles.  The  present 

model  is  simply  a  model  by  which  the  procedure  of  manipulations 

typical  for  our  schemes  of  approach  is  demonstrated.  It  is 

too  hasty  to  assume  that  the  model  simulates  either  the  ion 

system  or  the  electron  system  in  a  real  plasma.  As  is 

discussed  in  subsection  4.  10,  it  is  necessary  to  investigate 

model  II  as  is  given  in  the  next  section  for  evaluating  more 

precisely  all  the  effects  of  possible  interactions  on  the  evolution 

4 

of  the  ion  distribution  as  well  as  of  the  electron  distribution. 

A  te3t  particle  may  experience  three  main  sorts  of  interaction 
with  the  field  particles.  1)  Strong  interaction  with  its  nearest 
neighbor;  2)  Interaction  with  particles  which  are  not  nearest 
neighbors  but  are  fairly  _j.ose  to  it;  each  of  those  field  particles 
has  its  own  nearest  neighbor  and  exerts  a  fluctuating  force  to  the 
test  particle.  3)  The  force  exerted  by  remote  field  particles; 
the  force  due  to  each  field  particle  is  trival  but  culminates 
to  a  Vlasov  type  force  in  the  macroscopic  sense. 
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Let  us  take  particle  i  for  the  test  particle;  by  taking 
2  for  v  f  Eq.  (2.  4)  yields 


U) 


V3  j/k 


{<?ij  +<5?ik) 


a 

•sir 


*(3) 


(ijk)  dXjdXj^  = 


=  0 


and  Eq.  (2.  6)  leads  to 


(4.1) 


t)  ty-  /v. 


&'ik] 


P(3)(ijk) 


(4.2) 


Here 

F(3,(ijk)  =  V3  j  D(N)TTd^  (l^i,j,k) 

=  F(3)(X.X.Xk;t)l  if  i  <j  <k 
=  F^X.X.X^  t),if  j  <i  <  k 

etc. 
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Equation  (4.  1)  i.  a  partial  differential  equation  of  F  1  <V> 
with  X(  and  t  for  the  independent  variables .  It  is  possible, 
in  general,  to  determine  F(1)(X.;t)  along  a  certain  trajectory 
during  the  time  period  from  t  to  t  +  s.  If  the  trajectory  is 


specified  by 


.  t  i  f 

X.  =  X.  +  <s),  X.  =  pif 


(4.4) 


x.  (0)  =  0 

l  • 


p.  (0)  =  Pi.  r.(0)  =  r. 


the  function  is  given,  along  the  trajectory,  by 


F(1)(x|  ;  t  +  s) 

Considering  the  above,  we  coarse-grain  Eq.  (4.  1)  with 
to  s  from  0  to  t  along  trajectories  determined  by 


I 

P, 


(independent  of  s) 


t 


dt  m 


(4.5) 

respect 


(4.6) 


rl  <°>  =  ri 


(trajectory  0  ) 


From  now  on* trajectories  determined  by  Eq.  (4.  6)  will  be 
referred  to  by  class  0,  and  we  assume  that  t  is  much  shorter 
than  the  macroscopic  time  scale  and  much  longer  than  the  time 
scale  of  interaction  between  two  neighboring  particles.  The 
coarse -graining  of  the  first  two  meml'  rs  of  Eq.  (4.  1)  leads  to 


t 


3  s  m 


— F(1)(X.\  t  +  s)  ]  ds 
^ri  along  0 


=  —  J  _  ur  F  (-V +  s)ds 


(4.7) 


where 


<  F(1,(i)>  =  ~  f  F^fx'jt  +  s)  ds 

0  T  J  o  1 

(Along  trajectories  0)  (4.  8) 

The  integral  term  in  Eq.  (4.  1),  by  the  same  coarse-graining, 
yields  7* 

i  r  , 


Jo  j 


(4.9) 


Note  that  the  time  average  is  made  after  the  summation  2...  • 

JK 
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Here  X'  .  is  a  function  of  3  according  to  Eq.  (4.  6}  and 

1  JK 

means  that  similar  f  ctors  are  to  be  summed  up  with  respect 
to  (N-l)/2  pairs  of  particles.  See  (2.  5).  It  is  noted  that 
the  factor  of  a  pair  cannot  be  decomposed  into  two  factors, 
each  made  with  respect  to  one  single  field  particle,  after 
the  coar se-graining  with  respect  to  time;  the  reason  is  that 
the  interaction  among  three  particles  i,  j  and  k  is  a  non-linear 
phenomenon.  Hence  the  equation  of  evolution  of  <  F^(i)  >  q 
is  given 


•SF7'< 

1 


F(1)(i)> 


+  <  J>  =  0 
0 


(4.  10) 

In  order  to  obtain  F  ^  to  be  inserted  .1  the  above 
integrand,  it  is  necessary  to  solve  Eq.  (4.  2),  a  three  body 
problem.  The  solution  must  be  given  along  trajectories  0 
determined  by  Eq.  (4.6).  In  the  following  howe’  er,  we  obtain 
particular  solutions  of  Eq.  (4.2)  integrated  along  several  different 
classes  of  trajectories  of  particle  i,  j  and  k  by  successive 
approximation  methods.  Later  the  differences  among  the  results 
due  to  the  differences  among  these  classes  of  trajectories  are 
eliminated  by  consideration  of  proper  coarse-graming  operations 
with  respect  to  similar  particles.  For  the  purpose  of  providing 
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during  the  tim?  between  t  and  t  +  s  we  rewrite  £q.  (4.  2) 
as  follows: 


>  K  * 
*'i  *  ■  *'i  * 


♦ 


A*  ^ 


<->- 1  1  £  o* 1  ry  *  (z) 

<J}i  'Ji *>•—  *<Ai  |F  (![iYi|t**>  *° 

0  P j  d  ?k  J 


/V  /  ^ 


(4.  11) 


As  is  stated  in  Appendix  A,  is  invariant,  if  Eq.  (4.  11) 

is  integrated  along  a  trajectory  specified  by  XV  ,  XV  ,  X'^ 
obtained  as  functions  of  s  by  solving  Eqs.  (A.  1).  Since 
solution  of  Eqs.  (A.  ’)  is  a  three -body  problem,  we  consider 
combinations  of  two-body  problems  which  approximate  the 
three-body  problem  under  various  circumstances.  For 
each  case  of  approximation,  we  take  a  different  set  of 
differential  equations  in  place  of  Eqf.  (A.  1).  For  each 
case,  the  inifial  conditions  are  always  given  by 


x'.  =  X  ,  x'  =  X  ,  x'  =  X  (4.  12) 

1  J  j  K  JC 

8=0 

Each  different  set  of  differential  equations  replaced  for  Eqs.  (A.  1) 
determines  a  differei  ^et  of  trajectories.  Along  those 
trajectories,  F^  is  no  more  invariant;  it  is  a  function  of  s. 

In  each  case,  F^  is  obtained  by  solving  Eq.-  (4.  11)  by 
a  different  sucessive  approximation  method.  In  order  to 
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introduce  proper  sets  of  differential  equations  in  place  of 


Eqs.  (A.  1),  we  define,  as  is  detailed  in  appendix  B,  W„ 
which  is  the  probability  that  particle  i  and  particle  j  are 
their  mutual  nearest  neighbors.  Similarly  we  may  define 


and  W. 


jk 


WijWik  a  °.  3  /  K  '  ' 

The’-.,  defined  by 

F(i/j/k)  =  (1-W..)  (1-W.k)(l-W.k)  F{3)(ijk) 

*  -  wik  -  vF<3><ijk| 

(4.  14) 

F(i/j/k)  is  the  joint  distribution  of  particle  i,  particle  j  and 
particle  k  among  which  none  is  probably  the  nearest  neighbor 
of  the  others.  Similarly  we  define 
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y (ij/k)  =  w  F(3)(ijk) 

F(ik/j)  =  W.,  F(i)(ijk)  (4,  15) 

lk 

F(i/jk)  =  W  F(3)(ijk) 

JK 

Obviously 


F(i)(ijk)  =  Fi/j/k)  +  F(ij/k)  +  F(ik/j)  +  F(i/jk) 

(4.  16) 

In  view  of  the  definition  of  the  W' s,  those  members  in  the 
right  hand  side  of  Eq.  (4.l6)cannot  be  significant  all  at  the 
same  time;  If  one  is  significant,  the  others  are  trivial. 

In  solving  Eq.  (4.  11)  for  obtaining  F^(ijk)  appearing  in 
(4.  14)  and  (4.  I5),  we  have  four  different  methods  of  approximation. 
For  example,  for  calculating  F^(ijk)  involved  in  F(i/j/k),  we  may 
assume  that  interactions  among  the  three  particles  are  always  weak, 
ar.d  hence,  in  the  zeroth  approximation,  those  three  particles  are 
regarded  as  in  free  flight.  On  the  other  hand,  for  F^(ijk)  in  F(ij/k), 
it  is  essential  to  assume  that  the  interaction  between  particle  i  and 
particle  j  is  strong,  although  these  two  particles  interact  with 
particle  k  only  weakly.  On  obtaining  F^(ijk)  in  four  different 
cases  respectively  by  four  different  methods,  the  results  are 
presented  as  follows: 
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i-  r(3)(ijk)  =  F(i/j/k)(0)  +  F(i/j/k)(1)  +  -  -  - 

II.  F(3)(ijk)  =  F(ij/k)(0)  +  F(ij/k)( 

II*  F(3)(ijk)  =F(ik/j)(0)  +  F(ik/j)(1)+  -  -  - 

III.  F(3)(ijk)  =  F(i/jk)(0)  +  F(i/3k)(1)+  -  -  - 

After  obtaining  F^(ijk)  in  four  different  cases, 
we  substitute  each  of  the  solutions  to  its  proper  place 
in  (4.  14)  and  (4.  15).  Thus  for  F^(ijk)  as  presented 
by  (4.  6),  we  obtain 


F(3){ijk)(0)  +  F(3)(ijk)(1) 


+. 


where 


F(3)(ijk)(0)  =  (1  -  W..  -  W.k-  W.k)  F(i/j/k)(0) 
+  W..  F(ij/k)(0)  +W  F(ik/j)!0) 

1 J  IK 

+  wjk  Fa/3k)<0> 


(4.  17) 


{4.  18) 
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F(3)(ijk)(1)  =  (1  -  W  -  V.k  -  W  )F(i/j/k)(1) 


+  W..F(ij/k)(1)  +  W.kF(ik/j) 


(4.  19) 


(1) 


+  W.kF(i/jk)(1) 


etc.  (4. 20) 


In  the  next  step  of  our  treatment,  we  substitute 
F^(ijk)  presented  by  (4.  18)  in  the  integral  member  in  Eq.  (4.  1) 


J  =  jk  JJ  **>*** 


(4.21) 


where  J^,  for  example,  is  due  to  F^(ijk)^  given  by  (4.  18): 


r(°)  ,  y(0)  .(0)  ^  .(0) 

J!  +  J2  *  J2 1  +  J3 


(4.  22) 
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(4.23) 


<£3  -^ki-jr;  «1^(iV3>(0,"j«k 


(4.24) 


(4.25) 


Similarly  we  have 


J 


(D 


(4.27) 


in  which 


(1)  = 


1 

V2 


jk 


Jj,^ 


+  ^ik> 


a 

5  Pi 


(  1 


w..  .  w.u 

IJ  ik 


Wjk )  Fti/j/k)*1*  dX.dXk 


(4.  28) 


J2“’  =  7JjJ  «*,  ♦?*>  ■  -sir  wijF<i)/k)a>“j“k 

^  (4.29) 


Jz-I,=-£-£  II  wikF<ik/j)<1>«jdxl 


(4.30) 
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J 


3 


(1) 


1 

=  -v? 


(4.31) 


etc. 

For  coarse-graining  J  with  respect  to  time  from  t  to  t  +  t. 

we  may  have  more  than  one  method.  In  (4.  9)i  X'  ^  changes  as  a 

function  of  s  according  to  (4.  6),  while  X'  .  and  X'.  are  invariant 

J  K 

during  the  time  of  coarse-graining.  If  we  choose  for  X' .,  A1  .  and  X' ^ 
a  different  class  of  fui.  •  ms  of  s,  we  may  obtain  a  different  result 
of  coarse-graining.  In  the  following  we  obtain  four  different 
functions  for  F^(X'.  X'  .  X' t  +  s)  by  integrating  Eq.  (4.  11)  along 

i  j  K 

four  different  classes  of  trajectories.  These  classes  of  trajectories 
are  also  different  from  the  one  determined  by  Eq.  (4.  6).  The  result 
of  coarse  graining  J,  and/or  a  part  of  J,  along  a  class  of  trajectories 
is  different  from  the  result  made  along  another  class  of  trajectories. 
Later,  however,  it  will  be  shown  that,  by  a  proper  coarse-graining 
operation  with  respect  to  similar  particles,  the  differences  among 
the  results  due  to  the  difference  among  classes  of  trajectories  along 
which  the  coarse -graining s  with  respect  to  time  are  done,  are 
eliminated.  In  the  following  we  obtain  for  F^  four  different  functions 
of  s,  respectively  corresponding  to  four  different  ser.c-s  expansions 
given  by  (4.  17). 
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u .2.  Equations  for  sub-systems. 

4.2.1.F(3)(ijk)  in  F(i/j/k) 

According  to  (4.  17),  we  define  F(i'  / j '  /k'  ;t  +  s)' 
for  F(i/j/k  )^  where  X.,  X,  XR  and  t  are  replaced  respectively 

by  X!  ,  X'.,  X',  and  *  +  s.  Here  X’.,  X'.  and  X'  are  functions  of  s 

1  J  (0) 

as  will  be  specified  soon.  We  assume  that  F(i‘  /j'  /k'  ;  t  +  s) 

satisfies 


iyx) 


P; 


^ds  +  m  '  ?>t. 


+1L.JL,  +  _JL-  .JL,)  F(i7j'/k';  t  +  s){0)  =  0 


m  r\T  . 

J 


(4.32) 


and  F(i’  /j'  /k'  ;  t  +  s)^  for  u  >  0  satisfies 
P.  /  „  P:/ 


Pk'  a 


(A.  +  _J_1  .  +  AL  .  -2_  +  —  .-2-r  )  Ffi'/jfk' ;  t  +  s) 

'  d8  m  dr.'  m  dr-  m  dr  k 


(u) 


=  ■  [<31/  +J'ik,)-a|T:  *  ^i#  +  dV>  ‘sp( 

“  J 

+  (?ki  +v?kj  )  •  sfc  1  F(i'/j7k';  1  +  s)^  l) 


u  =  1.2.  3.--- 


where^.'  is  a  function  of  r.'  -  r.'  instead  of  -  r..  By 

^  J  J 


taking 


F(i'/j'/k'  ;  t  +  s)(0)  =  F(i'  ;  t  +  s)(0)F(j' ;  t  +  s){0) 

F(k';  t  +  s)<°> 


(4.33) 


(4.  34) 
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equation  {4.  32)  leads  to 


~j)  F(i';  t  +  8)(01  =  0 


FU'it  + 

TTk>  r<k  •  ‘  +  5 


(4.35) 


On  integration  of  Eqs.  (4.  3)  along  trajectories  determined  by 


dr  / 
.1 


drk' 

~ai~ 


/ 


m 


Pi 


/ 


=  P; 


P 


k 


/ 


=  P 


k 


/ 

r. 

i 


(0) 


r  , 

i 


r.'(0)  = 


r .  , 
J 


rk,!))  ■  rk 


[  trajectory  I] 


(4.36) 


(4.37) 


we  have 
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F(i'  ;  t  +  s)<°>  =  F(‘.,  t) 
F(j'  ;  t  +  s)(0)  =  F(j  ;  t) 
F(k' ;  t  +  s)'0)  =  F(k  ;  t) 


(0) 

(O' 

(0) 


(4.  38) 


The  next  procedure  is  to  obtain  F(i'  /j'/k1  ;  t  +  s)^  by 
solving  Eq.  (4.  33)  for  4=1.  On  substitution  of  (4.  38)  in 
the  right  hand  side  of  Eq.  (4.  33),  we  have 


<TT+4r  -sir-  +4--  TFT'  +4NF(i'/j'/k':t  +  8) 


"  t  ^ij  +<3*ik  }  •■^7'  +  ^ji  +i>jk) 


(1) 


+  fa  +^j'>-  1  F(i'/j7k';,+  s)<°) 


(4.  39) 

Since  the  right  hand  side  member  is  already  given,  we  may  obtain 
F(i*  /}'  /k*  ;  t  +  s)(1>  along  trajectories  .  I  determined  by  Eqs. 

(4.  36)  and  (4.  37).  We  may  repeat  the  same  procedure  to 
obtain  F(i’  /j1  /k1  ;  t  +  s)^  for  n  >  1.  It  is  obvious  that 


F(i'  /j'  /k'  ;  t  +  S){U)  =  0, 
if  u  i  0,  s  =  0 


(4.40) 


31 


since  w«  assume  that 

F(i/j''k;t)(0^  =  F'!)  (X.  ;t)  F{1)(X.;t)  F(1)(X,  ;t)  (4.41) 

1  J  K 

4.2.2.  F(i)(ijk)  in  F(ij/k) 


In  this  cae,  particle  i  and  particle  j  are  in  a  close 

interaction.  By  taking  the  second  expansion  of  F^(ijk)  given 

(0) 

by  (4.  17)  and  defining  F(i‘  j1  /k1  ;  t  +  s)  ,  etc.  in  the  same  way 
as  of  F(i'  /j'  /k'  ;  t  +  s)^,  etc.  we  have 


( JL_  +  _!/_  .  JL.  +  Jj_  .  JL_  +  .  _JL_ 

ds  m  dr./  m  dr./  m  dr,  , 
1  1  ^ 


•t+»>(0)=0 


rJ_  +  JV_  *  ■  pi'  £_  +Zlll  J_ 

Lds  m  '  dr/  m  ’  dr'  m  '  d  r,  ' 
•  J  K 


~r  /  >  9  *  \ 

+  ,Aj  •  "5p~’  +/i  -4:0  F(iV,k';t+3) 


(4.42) 


-  .  r  %  3 
1  ^ik  -Tp, 


+  3>'  b 

J*  ’  *PT' 


32 


.  (U-U 

A_,]F(i*j'/k';  t+  s) 

°pk 


4  =  1,2,3,... 

Those  equations  are  integrated  successively  along  trajectories 


(4.43) 


II  determined  by 


(4.41) 


33 


•IlillillllllMItl' 


We  may  assume  that 


F<i'j7k'  :  t  +  s)(0)  =  F(i'j':  t  +  s)(0)  Fik'.t  +  s){0) 

(4.46) 

where  F(i'  j'  ;t  +  s)^  and  F{k'  ;t  +  sf  are  solutions  of 


#JL  +  ll.  .  +Zil  5  , 

'gs  m  ar.  m  a** . 

i  3 


•  6  + 

*  "ap“^  r  "ji  oPj 


3**  •  ^  )  F<i'i/;  1  r  s)(0)  “  0  (4-47) 


and 

a  Pk'  <0' 

(It  +  -sr  -rl-0  F(k'; 1  +  s)  =  0  <4-48> 

ds  m 

when  integrated  along  rajectories  II  determined  by  (4.40) 
and  (4.45).  Since  we  assume  that 


when 


F(i'j7k' ;  t  +  s)(u*  =  0 
|i  >  0 
s  =  0 


it  is  easily  seen  that 


when 


F(i’  j'  ;  t  +  s)<°>  =  F(-l)(ij) 
F(k' ;T  +  s)(0)  =  F(1)(k) 

s  =  0 


(4.50) 


(4.51) 
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Solution  of  ’Sq.  (4.47)  is  a  two  body  problem;  we  may  imitate 

the  method  of  Boltzmann  setting  Stos szahlansatz  for  the  initial 

condition  of  a  binary  *  ollision.  Of  course,  it  is  to  be  expected 

that  the  solution  will  become  meaningless  as  the  distance  between 

iwo  particles  i  and  j  increases.  However,  at  the  same  time 

W.  .  tends  to  vanish,  and  hence  such  meaningless  solutions  do 
ij 

not  affect  the  final  result,  F(ij/k). 


4.0.-.  F(i)(ijk)  in  F(ik/j) 


We  may  repeat  the  same  treatment  in  this  case  as  in  the 
case  of-  F(ij/k),  simply  by  changing  subscripts.  The  trajectories 
which  correspond  trajectories  II  are  denoted  by  trajectories  II1 
in  this  case. 

4.?.^.F(i)(ijk)  in  F(i/jk). 


Trajectories  along  which  coar se -graining  operations  are  done 

are  named  trajectories  III.  The  treatment  is  similar  to  those 

in  the  two  preceding  cases  under  the  condition  that  the  interaction 

oetween  j  and  k  are  strong. 

Considering  F^(X'  .  X'  ^X'  t  +  s)  obtained  in  fouv 

cases  in  the  above,  we  may  coarse  -  grain  J  J.^,  J,.  ^ 

(0)  ILL 

and  J,  ,  with  respect  to  time,  respectively  along  trajectories 


35 


I,  II,  II'  and  m. 


4.  3.  Coarse-graining 

Coarse -Graining  Jj^,  with  Respect  to  Time. 

According  to  (4.  22),  (4.  33)  and  (4.  36),  we  have 

<j1<0)>1  ■  e  -f-  r0T  r  J J  <£/+#k>  w:  <*-  w'ii  -  WV  V 

jk 


F(i';  t+  s)(0)F<j\  t  +  s)(0)F(k/;t+  s)(0)dX'j  dX^] 


£  rI 1 1  -  wij -  wik-  V  (^J  +^k>  F<  l,«)  F"’<k> 


-gr~  F<1)<i)  dxj  dxk 


I“  wij>Jlj  f<\)  dxr4_  r(I’(i) 


(4.  52) 

In  the  above,  it  is  noted  that  £  means  the  summation  with 

jk 

respect  to  a  possible  set  of  pairs  made  of  all  the  particles 
except  for  the  test  particle;  hence  the  number  of  the  pairs  is 
(N-l)/2.  By  the  time  <  verage,  the  symmetry  is  produced  in 
the  distribution  function.  It  is  also  assumed  that 
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va  ik 


ff  W-.  F(1)(j)  F(1)F(J)(k)  dX  dX. 

<j  j  ik  J  K 


i.>  negligible  as  compared  with  the  other  members.  <  Jj  >  gives 
ine  Vlasov  effects.  With  necessary  cautions  similar  to  those 
lor  (4.  5 L),  we  obtain 


1  N  - 1 


—  <  fjT.  V . .  .  — — F(Z)(ij)  dX. 

I  ’J  J  II 


II  V 


(4.53) 


<v 

c  "ir  v  l 


1<  IK  A  FtZ)(ik)  dX.^TT 

^  J  b'l}.  np.  '  '  k  II 


+  — - j  1  )(j)dX.  .  — —  j  F{Z)(ik) 


wikdXk': 


(4.  54) 


it  is  noted  that  j*  F(1)(X)dX/V  1,  and  /  JP  /- .v 1$^  m  (4*  5i) 

and  in  (4.  54).  Therefore,  the  second  member  in 

each  is  ignored,  and  the  sum  of  them  leads  to 


T  (0)  ,  T  (0) 

Ji  >11  f 


--  <  fX  W,  •  FU)(ij)  dX.-. 


•  U  U  *P. 


(4.  55) 
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This  will  be  converted  to  the  Boltzmann  collision  integral  by 
proper  treatment.  Here  the  convergency  of  the  result  is  assured 
by  the  existence  of  W. ..  Finally  we  have 


<  J 


3 


(  N-l) 

— z —  < 


(2) 

F  (jk)dX  dX  > 

J  k  III 


d  -(1),.. 

ftp-  <  F  (l)  > 

BPi  III 


(4.  56) 


It  appears  that  the  effect  is  similar  to 
involves  W.,  .  Therefore 


The  difference  is 


that  <J. 


jk- 


(0) 


(0) 


<3  *  >  <  <  J,  > 

J  III  1  I 


(4.  57) 


Hence  we  n  ay  ignore  the  effect  given  by  (4.  56) 


*.3.2.  Coarse -Graining  wit^  Respect  to  Time 

The  J^'  s  given  by  (4.  28),  (4.  29),  (4.  30)  and  (4.  3  1)  are 
coarse-grained  in  the  same  way  as  the  J^'  s.  The  trajectories 
along  which  those  functions  are  coarse-grained  are  respectively  the 
same  as  those  along  which  the  J^'  s  are  coarse-grained. 

<  There  are  correlations  among  three  particles 

i,  j  and  k  in  this  approximation.  It  is  easily  seen,  however, 
that  those  particles  are  m  itually  remote.  According  to  (4.  39), 
we  obtain 

F(i / j *  /k'  ;  t  +  a){1) 
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+ 


F(i'/j'/k'  ;  t  +  s)(0)da 


(4.58) 

Since 

F(i'  /j*  /k'  ;  t  +  s)(0)  =  F(1)(i;t)  F(1)(j;t)  F(1)(k;t) 
and  those  trajectories,  along  which  the  coarse-graining 
is  done,  are  free  m.ituall>,  the  summation  r  ^  results  in 
cancellation  among  the  effects  of  those  particles.  Hence  we 
ignore  the  present  effect. 


<r  ^  ^  +  c  J^,,  #  It  is  easily  shown  that 

integrations  of  Eq.  (4.43)  £  .id  me  similar  equation  for 
F^J  (i *  j '  /k*  ;  t  +  s)  along  trajectories  II  and  II'  lead  to 


J 

L  II 


(i).  _  n  - 1  ifrrrr'£/  a 

z-  ir  ‘  LJJ</ik'  W 


s  I 


Wij'  r  7ik  *  *fr  F<2)(i,j,;t+s)F(1)(k‘;t+s)d8dX'dX,Jd8 


"ij  ~  r(2)<i’J'it«) 


N-  1  _1_ 
~2V2 


Jo 


F(1)(k»;t+s)dedXJdxJ 


ds 


(4.59) 
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In  the  above,  the  first  member  in  the  right  hand  side  gives 
the  diffusion  effect  in  the  momentum  space  of  particles  i  and 
the  second  member  the  friction.  The  significant  feature  of 
those  terms  is  that  the  correlation  period  is  regu  lated  by 
W^j,  the  probability  of  strong  interaction  between  i  and  j. 


<  J 


(D^ 

3  III 


.By  means  of  similar  method  as  before,  we 


obtain 


<  J  (»-.  «ll 

III  2VS 


-g|—  ;  t+s)F(1)(i*;t+s)ds  dX'dX^Jds 


2V 


N-J  J__  r  r  if  7  '  _L_  W  ' 
T  J  nl  JJ  ^  ij  ^p7  jk 


S  'J’l'i  'W- t  F<2l(jk,;  ds 

0  J 


(4.  60) 


In  effect,  there  is  no  difference  between  <J?^  ^  +  'J 

^  ii  2  II 1 

(i) 

1  in' 


(i). 


and  <r  J . 
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In  summation  of  all  those  effects,  we  may  write 


<J 


U>.  =  .  (N-1)  j _  r  r  rr^t  '  ._i_  w  ' 

v3  T  J  ,p.' 


1J 


+  J  8  Jik  --^T  r(2)(i' j'  ;  t  +  8)  F(1)(k-  ;  t  +  s)  a 


8 


dX'  .  clX'  .  i  da 
J  k 


(N-1) 

V® 


h 

a  Pi 


F(Z)(i'  j'  ;  t  +  s)  F(1)(k'  ;  t  +  s)ds 

dX.'  dX.  '  I  d» 

J  k  J 

(4.61) 

where  it  is  noted  that  we  have  ignored  the  difference  among 
trajectories  over  which  those  coarse-graining  operations 
are  considered.  As  is  shown  below,  however,  the  difference 
is  not  essential  in  the  result.  The  effect  given  by  (4.61) 
is  of  the  Fokker -Planck  type. 
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*♦•3.3,  Distribution  functions  in  generic  spaces 


So  far,  various  functions  constituting  J  have 
been  coarse-grained  along  different  classes  of  trajectories. 

For  completing  Eq.  (4.  10),  however,  J  is  to  be  coarse-grained 
with  respect  to  the  particular  class  of  trajectories  named  0  as  is 
given  by  (4.  V,  It  is  necessary  to  show  that  the  differences  among 
those  classes  of  trajectories  along  which  J  and/or  parts  of  J 
have  been  coarse-grained  do  not  result  in  any  inconsistency 
in  the  final  result.  For  this  purpose,  we  compare,  for  example, 
the  results  of  coarse -graining  F^(X.  ;t)  along  different  classes 
of  trajectories  in  the  following:  According  to  (2.3), 


F(1)(X.;t)  =  VI*  D<N>  tt  dX. 

1  J  •  /  .  1 

j  n  J 

=  V6(X.-  X.*(t))  (4.62) 

* 

Here  X^  (t)  is  a  complicated  function  of  time  because  of  the  interaction  of 
particle  i  with  the  other  N-l  particles. 

First,  let  us  coarse-grain  F  ^  along  trajectories  determined 
by  (4.  6)  over  a  time  period  from  t  to  t  +  t: 


<  F  (1)(X  ;th  =  -1 

11  T 


J  F  (1)(X.'  ;  t  +  s)ds 


(4.63) 
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where 


X'.  =  X.  +  x.'(s) 
1  1  1 


Those  trajectories,  according  to  (4.  6),  are  independent  of  the 

other  particles.  As  is  shown  in  Fig.  4.  1,  none  of  those  trajectories 

is  similar  to  X.  '  (t  -t  3)  in  view  of  time  scale  t.  Hence 

<F  <1)(Xl;t)>  is  spread  over  a  comparatively  broad  domain  of  X^. 

There  may  be  a  group  of  v.  particles,  i  ,  i . . of  which 

1  1  c 


(1) 


,<H 


<FU  (X;t)>  =  <FJ^X;t)>  =  --- 
These  are  localized  in  a  narrow  domain  of  X.  We  define 
<F(1){X;t)>  by 

N  (1) 

N<  F(1)(X;  t)>  1  <  F .  (X; t ) > 

i  =1 


(1) 


=  <  F.,  (X;t)-> 


♦  v."  <F{11  (X;t)^  +  — 


(4.  64) 


j'  +  "i«  +‘  '  •  a  H 


(4.65) 
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Unlike  <F^^(X;t)%  ,  <F^(X;t)>  is  spread  mere  uniformly  over 
the  entire  domain  of  X. 

We  may  also  coarse-grain  F.^(X^;t)  along  another  set 
trajectories 


X'!  =  X.  +  x'!  (s) 

t  1  t  ' 


which,  for  example,  is  determined  by  (  4.4%)  % 


«  F. 


(1) 


(X.-t) 

i 


t  +s)  dS 


(4.66) 

Those  trajectories,  according  to  (4.4-l^are  determined  by 
taking  into  account  the  nearest  neighbors'  interactions. 

Therefore,  a  frajecotry  coincides  with  X.  (  t  +  s)  in  the  time 
scale  t.  Hence  «  F^^(X.;t-v>  is  much  more  localized  than 
<F,^,(X^',t>  as  is  illustrated  in  Fig.  4.  1. 


,X.; 


4  «F.(1)(X.; 
1  1 


t  )  » 


Nevertheless,  it  is  obvious  that 


(4.6  7) 


F 


(1) 


(X.  ;t)>  dX,  = 
i  i 


«F.(1)(X.;t)»dX. 


(4.  6h> 
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In  the  '•bove,  we  have  seen  coarse-graining  operations 
in  the  sense  of  (4.64)  and  (4.69)  -re  essential  for  reconciliation 
(4.  70)  betv'een  two  classes  of  time  averages.  In  other  words, 
time  averages  are  not  enough  for  deriving  proper  kinetic  equations; 
it  is  necessary  to  make  averages  over  many  particles. 

Applving  the  same  consideration  as  above,  we  make  average 
of  Eq.  (4.  10)  with  respect  to  particles  so  that  <  J>*  s  made  with 
respect  to  different  classes  of  trajectories  are  equivalent  to  those 
made  with  respect  to  trajectories  0.  By  this  procedure  of 
coarse-graining  over  particles,  we  obtain  in  effect  the  same 
function  as  f  defined  by  (Z.  7).  ^The  Boltzmann  equation  is  the 
very  equation  which  governs  the  evolution  sf  f  of  a  gas  where 
the  collisions  are  assumed  to  be  binary.  )  After  the  present 
coarse-graining  operation  with  respect  to  similar  particles, 
the  erfect  given  by  (4.  55)  may  amount  to  the  Boltzmann  collision 


integral;  because  of  W.^  in  the  integrand,  the  effective  interaction 
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terminates  quickly  as  the  dist.  c  between  particle  i  and  paiticle  j 
increases. 

4  Evaluation  of  the  Interaction  Effects. 

Because  of  the  W  s  involved  in  thost  interaction  integrals,  the 
integrals  converge  always.  As  is  given  in  appendix  A,  we  take  for  W 


W  =  exp  (  -9r  n  r^/4) 

where  n  is  the  local  number  density  and  r  is  the  distance  between 
the  test  particle  and  a  f!  ild  particle  under  consideration.  The 
range  of  space  of  a  close  correlation  (almost  binary)  ,  s  of  the 
order  of 

-1/3 


and  the  time  scale  is 


-  =  Rn/('p^/m) 

Here  <p>  is  the  average  magnitude  of  momentum  of  a  particle 


1 


If  r  R  , 


_e_e  =  o  riSli  7 

~n  -  m 

then  W  i.  The  Boltzmann  type  interaction  effect 
J ^ i  ^  is  of  the  order  of 
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0[  <J 


.(1), 


n  m 


-  TT 


n  R  ?  (  3kT  )TF(1)(i) 


n  '  m 


The  order  of  the  interaction  effect 
is  sh  wn  to  be 


oo 


nr  T  (!)  .  T  (D-,  Rn  r  4-e4r»  .  1  ^{1  ... 

Of  J/  >  +  ']  = - k  n  - - -  ar-— F'  '(x) 

1  1  <3kT/m)  J  Rn  r  P 


4- 


nc4 


3  mr  (3kT/m>ys> 


F(1)(i) 


Hence  we  obtain 


°r  j2(1)-+ -j2.(1)>i 

ot< a0[,-^r 


>a] 


(4.  71) 


4  19 

If  we  take  T  =  10  ,  n  =  10  ,  we  obtain 


ea  n 


kT 


-)  =  10 


-1 


( 


e 


kT 


) 


decreases  further  as  T  increases  and  n  decreases. 
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<  J3^>  is  the  flame  order  as  of  +  <  ^2'  • 

The  ratio  given  above  suggests  that  the  effect  or  remote 

interactions  may  be  ignored  in  ordinary  cases  where 
4  19 

T>  10  and  n  <10  ,  This  conclusion  appears  to  be 

different  from  those  obtained  by  other  authors,  based 
on  the  original  interpretation  of  the  BBGKY  hierarchy. 

It  is  noted  that  the  present  definition  of  close  interaction 
according  to  the  nearest  neighborship  is  different  from 
the  usual  definition  of  close  interaction  according  to 
the  polarization  effect  of  the  Debye-Huckel  type. 

One,  who  has  been  familiar  with  theories  based 
on  the  BBGKY  hierarchy,  might  ask:  Why  is  the  effect 
of  field  particles  which  are  not  nearest  neighbors  of  the 
test  particle  so  small  by  the  present  theory?  The  answer 
is  as  follows:  By  a  theory  based  on  the  BBGKY  hierarchy, 
the  test  particle  and  a  field  particle  may  interact  with  an 
indefinitely  long  correlation  time,  in  spite  of  the  presence 
of  perturbations  by  third  particles.  On  the  other  hand, 
the  behavior  of  a  particle,  by  the  present  theory,  appears 
in  microscopic  order  only  for  a  short  period  of  time  during 
its  interaction  with  its  nearest  neighboring  particles.  In 
other  words,  the  time  -<cale9  of  microscopic  orders  are  always 
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finite.  Note  that  the  above  statement  is  feasible,  only  because 
laws  of  interaction  forces  are  not  linear  with  respect  to  inter¬ 
particle  distances. 


4.  5.  Validity  of  the  Model 

As  is  investigated  in  section  III,  the  time  scale  of  a 
close  interaction  between  two  ions  is  much  longer  than  the 
the  time  scale  of  a  close  interaction  where  an  electron 
participates.  Therefore  the  present  model  appear  to  be 
suitable  for  simulating,  in  an  approximation,  evolution  of  the 
ion  distribution  in  a  real  system,  if  the  number  of  the  ions 
is  the  same  as  the  number  of  the  electrons  (charge  of  an  ion 
being  the  same  as  of  a  proton).  This  is  beqause  the  electrons 
appear  to  make  an  almost  uniform  background  in  view  of  the 
motion  of  an  ion.  But  at  the  same  time,  we  find  no  definite 
reason  for  saying  that  the  effect  of  electron-ion  (close) 
interactions  on  the  evolution  of  the  ion  distribution  is  negligible: 
although  the  average  magnitude  of  the  momentum  of  an  electron 
is  much  smaller  than  the  average  magnitude  of  the  momentum 
ox  an  electron,  (the  ratio  being  about  1/50  with  respect  to  the 
lightest  ion),  the  frequency  of  electron-ion  collision  (close)  is 
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about  50  times  higher  than  the  frequency  of  close  ion-ion 
interactions.  If  we  wish  to  consider  the  situation  more 
precisely  and  reasonably,  both  ions  and  electrons  must 
be  considered  as  discrete  particles.  In  view  of  the  above, 
model  II  will  be  investigated  in  the  next  section. 
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V.  MOL^L  II 


As  stated  in  subsection  4.  10,  the  validity  of  model  I 
as  simulating  a  real  system  is  fairly  dubious.  It  may  be 
plausible  to  think  that  the  model  has  a  meaning  only  as 
providing  an  example  of  manipulation  of  reducing  the  Liouville 
equation  to  coarse-grained  kinetic  equations  governing  evolutions 
of  subsystems.  In  this  section,  we  consider  a  model  which  may  be 
more  realistic  than  the  previous  one.  The  plasma  which  the 
present  model  is  intended  to  simulate  is  the  same  as  is  considered 
in  the  last  section;  the  gas  is  fully  ionized  and  the  number  of  the 
ions  is  the  same  as  the  number  of  the  electrons. 

We  consider  four  particles,  two  ions  and  two 
electrons,  for  representing  the  field  particles.  By  this  way, 


it  is  necessary  to  treat  the  evolution  of  a  five -particle  distribution 
(5) 

function  F'  (X^X^X^X^X^jt)  where  one  of  the  five  particles  is  the 
test  particle.  However,  as  a  simpler  approach  we  shall  consider 
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where  particle  o  is  the  test  particle,  either  an  ion  or  an 
electron,  particles  i,  j,  k  are  electrons  which  represent 
the  N/3  groups  of  the  field  electrons,  and  particles  £  **  ,  r 
are  ions  representing  the  N/3  groups  of  the  field  ions;  2N 
is  the  total  number  of  the  field  particles,  N  electrons  and 
N  ions.  The  evolution  of  F^(o)  of  the  test  particle  is 
obtained  by  integrating  partly  the  Liouville  equation  of 
the  N+l  particle  system  as  follows: 


,  3  ,  *o  3  .  F 0),  . 

+—  *  Tr“)  F  (o) 
o  o 


+  ^7  ^  [Jj  «*,  ♦£>,>•£"  *<*><6. 


F(3|(o,k.,)dXkdX, 

ro 


+  •  S^F<3)(0,’1'''  1  dX„dXJ  1  =  0 

(5.2) 

After  coarse -graining,  r  -  will  be  replaced  by  multiplier 

ijk  *r\  r 

2N/6.  Equation  (5.  2)  corresponds  to  Eq.  (4.  1)  for  model  I.  Cor 
responding  to  £q.  (4.  2).  we  have  the  following  three  equations: 
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+  _p~_  J _ +  pg  a  h  pP  a 

*t  m  'Sr  m  ‘?r  m.  '  irft 
t  o  a  a  (>  p 


/>  0“ 


+  <4  +  +  <Jao  +  ^ap) 


rr  *> 


tft>o+?Ba>  '  -TST  1  r<3)«>’  “■  P>  =  0 


Po  vpa 


BPp 


(a,  P)  =  (i.  j).  (k,§  ),  (r,  ,  c  ) 


where 


m.  =  m.  =  m,  =  electronic  mass 

i  j  k 


m  =  m  =m  =  ionic  mass. 


The  interaction  to  the  zeroth  approximation  is  given  by 


(5.3) 


j(0,= ^  ,3k  [  j  J  ^  ^  '-k  f,3)<°^><0Vxj 


+  JJ<£k  F<3>(o,k,.)<°>dX1dXf 


+  IJ'5ro.t5oj'-^  F(3)(o,r.  )<%XJXr  1 
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=  J(0)(oee)  + 

\ 

J(D)(oee)„  +  J(0)(oee) 

s 

t  +  J^(oee) 

+  J(0)(oel)  + 

J(0)(oel)  +  J(0)(oel)  , 

+  J(0)(oel) 

p  P 

3 

f  j'0)(oll)  + 

t 

J(0)(oII)  +  J(0)(oII)  , 

3  3 

+  J(0)(oII) 

3 

r(0), 


(5.4) 


where  J'  '(oeI)^f  for  example,  is  an  effect  due  to  field 

parades  k  and  e  ,  k  representing  field  electron*  and  ?  representing 

field  •  ions.  Those  members  in  the  above  may  be  obtained  as 

analogous  to  those  J^1  s  obtained  in  the  last  section.  In 

coarse-graining  these  terms,  it  is  necessary  to  pay  attention 

to  the  difference  between  the  time  scale  of  ion-ion  interaction  and 

the  time  scale  of  those  in  which  electrons  are  involved.  The 

time  scale  of  ion-ion  interaction  is  denoted  by  while  the  time  scale  of 

those  involving  electrons  is  denoted  byje*  Then,  according 

to  the  investigation  given  in  section  III,  it  holds  that 


II 


If  the  macro acopis  time  scale  t  m  is  larger  than  t^,  we 


may  choose  j  so  that 


T  '■*  v>  T  >>  Ttt  ^>T 

m  lie 


(5.  5) 
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t  is  the  very  time  scale  with  which  we  may  coarse-grain 
those  terms  in  Eq.  (5.4).  All  the  procedures  of  coarse- 
graining  are  the  same  as  demonstrated  in  section  IV. 

Taking  an  ion  for  test  particle  o. 

<J(0,(Iee)j>  +  <J(0)(IeI)1>  +  <J(0)(III)1  s  (5.6) 

gives  the  effect  of  the  Vlasov  type  interaction.  The  effect  of 
the  Boltzmann  type  interaction  among  nearest  neighboring 
particles  is  presented  by 


<J(0)(Iee)3  >  + 

<J(0,(Iee)r,  '  > 

+ 

<J(0)(IeI)?  >  + 

<J(0,(IeI)„  ' 

+ 

<J(0)(III)3>  + 

<r  J(0)(III)a  '  > 

We  may  also  obtain 

<J^(Iee)2  >  V 

^J(1)(Iee)3'^ 

+ 

<J(1)(IeI)p  >  + 

<J{1)(lel)p 

+ 

<J(1)(III),  >  + 

<  J<1,(III)S  - 

(5.8) 
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for  the  significant  effect  in  the  approximation  of  the  next  order. 
By  comparing  the  effect  given  by  (5.  8)  with  the  effect  given  by 
(5.  7)  in  the  same  was  as  we  did  by  (**.  71),  we  may  see  that  effect 
(5.  8)  is  negligible  in  ordinary  cases.  The  situation  is  explained 
as  follows:  1)  First,  let  us  take  an  ion  for  the  test  particle  o. 

Then,  take  an  electron  for  the  particle  nearest  neighboring  to  the 
test  particle  o.  The  effect  of  the  interaction  (a  part  of  the 
zeroth  approximation  effect)  is  larger  than  the  effect  of  the 
force  exerted  by  two  mutually  neighboring  field  particles 
of  which  one  is  an  electron.  Note  that  the  time  scales  of 
two  events  are  the  same.  Take  an  ion  for  the  nearest 
neighboring  partide  of  the  test  particle  o.  Then  the  effect 
of  the  interaction  is  larger  than  the  effect  of  the  force  exerted 
by  two  mutually  neighboring  field  particles  of  which  both  are 
ions.  2)  The  situation  is  similar  for  the  test  particle  which 
is  an  electron. 
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"V.  MODEL  in 


Suppose  that  a  fully  ionized  gas  is  composed  of  N 
electrons  and  N'  ions  of  a  single  species.  Here  the  charge 
of  an  ion  is  denoted  by  Q,  and 

Q  =  v  e,  e:  electronic  charge 

N/N'  =  „  (4. 

If  x  is  much  larger  than  unity  and  we  wish  to  investigate  the 
evolution  of  the  distribution  of  ions,  the  models  given  so 
far  are  not  feasible:  Model  I  is  not  feasible  because  neither 
ions  nor  electrons  can  be  represented  b1  a  uniform  charge 
distribution  in  the  present  case;  model  II  is  not  feasible 
because  N  >>  N1 .  Due  to  fee  large  mass  and  the  large  charge 
of  an  ion,  two  nearest  neighboring  ions  which  are  separated 
at  a  distance  much  larger  than  the  average  distance  between 
two  nearest  neighboring  electrons  may  have  a  strong  mutual 
correlation  of  which  the  time  scale  is  much  longer  than  the 
time  scale  of  the  behavior  of  an  electron  interacting  with 
an  ion  or  another  electron.  See  Fig.  6.  1.  This  condition 
is  quite  favorable  for  applying  the  ordinary  BBGKY  hierarchy, 
up  to  the  third  equation,  to  the  investigation  of  the  evolution  of 
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the  ion  distribution.  Here  the  test  particle  is  an  ion,  the 
field  particles  appesring  in  the  first  equation  of  the  BBGKY 
hierarchy  are  also  ions,  and  those  field  particles  appearing 
in  the  integrand  of  the  second  equation  are  electrons.  The  test 
ion  and  a  field  ion  may  interact  through  a  space  in  which  many 
electrons  are  present  and  the  distribution  of  those  electrons 
are  polarized  by  the  field  induced  by  the  two  ions.  The  case 
has  been  investigated  by  many  authors^*  as  based  on  the  BBGKY 
hierarchy. 


VII.  CONCLUDING  REMARKS 


If  the  number  of  ious  is  the  same  as  the  number  of 
elec  -ons  ir.  a  plasma,  U.*re  are  two  basic  types  of  inter¬ 
particle  interactions  in  the  zeroth  approximation:  One  is  of  the 
Vlasov  type  and  the  other  of  the  Boltzmann  type  characterizing 
interactions  among  nearest  neighboring  particles.  In  higher  order 
approximations,  mutual  perturbations  among  those  basic 
interactions  resu;f  in  secondary  effects  which  are  not 
significant  under  ordinary  conditions  of  high  temperature 
and  of  low  density. 

As  the  ratio  betwv  -n  the  number  of  electrons  and 
the  number  of  ions  incre?,aes  in  a  real  system,  the 
simulating  model  becomes  more  complex;  in  cases  where 
the  ratio  is  extremely  large,  the  BBGKY  hierarchy  is  useful 
for  investigating  the  evolution  of  the  ion  distribution. 

The  presence  of  any  external  force  field  of  which  the 
time  scale  is  much  larger  than  the  maximum  time  scale 
of  close  inter-particle  interactions  does  not  affect  the 
validity  of  the  present  approach.  If  ti.^re  is  an  external 
force  field  of  which  the  time  scale  is  comparable  to  or  smaller 
than  the  time  scale  of  any  close  inter-particle  interactions,  the 
method  must  be  considerably  modified^. 
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APPENDIX  A 


Solution  cl  Eq.  (4.  11) 

Equation  (4.  11)  is  a  partial  differential  equation  where 

F^(X'  .X'  .X' ,  ;  t  +  s)  is  unknown  and  X1  ..  X'  X' ,  and  s 
i  J  *  1  J  k 

are  independent  variables.  It  is  possible  to  obtain  particular 
solutions  by  taking  arbitrary  functions  of  s  for  X'  ,  X1  ..  X'.  . 

s  J  1C 

If  one  takes  for  X'  X'  j  and  X'  ^  particularly  those  which  satisfy 
the  characteristic  equations  of  Eq.  (4.  1 1) 


m 


(A.  1) 


then  the  solution  of  Eq.  (4.  11), is  invariant  with  respect 

to  s.  Ii  is  noted  tha.  X1  .(s),  X'  ,(s)  and  X' ,  (s)  which  satisfy 
Eqs.  (A.  1)  present  a  set  of  trajectories  of  particles  i,  j  and  k 
which  are  in  mutual  interaction.  The  solution  of  Eqs.  (A.  1)  is 
a  three-body  problem.  In  order  to  avoid  the  difficulty  of  solving 
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the  three-body  problem  precisely,  one  may  choose  a  set  of 
approximate  differential  equations  in  place  of  (A.  1)  depending 
on  circumstmces,  and  may  integrate  Eq.  (4.  11)  along  a  set  of  tra¬ 
jectories  given  by  Xi  ,  X'^  and  X1  ^  which  satisfy  the  approximate 
differential  equations.  In  this  case,  the  solution  is  not 
precisely  invariant  with  respect  to  s. 


APPENDIX  B 

The  Probability  of  Two  Particles  being  their  Mutual  Nearest  Neighbors"1 

If  it  is  defini'e  that  n  particles  exist  in  a  space 
domain  of  unit  volume,  the  probability  of  no  particle  being 
present  in  a  domain  of  volume  v,  which  is  smaller  than  unit 
volume,  is  given  by 


wq(v)  =  exp<-nv)  i  B.  1) 

The  proof  is  as  follows:  Let  us  consider  a  variation  6v  of  v  , 
the  probable  number  of  particles  in  6v  is  n6v.  The  probab;'  „y 
of  any  v-  particles  appearing,  at  the  same  time,  in  6v  is  denoted 
by  w  (6v).  Then,  we  have 


n 

r  v  vv  (6v)  =  n6v  (3. 2) 

v 

V*  1 

Obviously 

n 

w(j(6v)  +  r  w  (6 v)  =  1  (B.  3) 

v  -  1 
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Since  we  know  that 
lim 

6v  0  w  (6v)  cc  (6v)v 
V 

(B.  2)  leads  to 

lim  w.  (6v)  =  n6v 
6v  0 

Hence  QB'  .  3)  yields 

w  (6v)  +  n6v  =  1 
o 

at  the  limit  8v  -*  0. 

On  the  other  hand,  we  have 

w  (v  +  6v)  =  w  (v)  w  (6v) 
o  o  o 

By  expanding  the  left  hand  side  of  (  B>  3),  we  obtain 
dw 

w0(v)  +— =  wo(v)  (  1  -  n6v) 

where  (  &,4)  is  taken  into  account,  or  we  have 

dw  (v) 

- 2 —  =  -  ndv 

w  (v) 
o'  ' 


(EL4) 


(B.  5) 
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On  integration  of  the  above  under  the  condition  wq{0'  =  1, 
we  obtain  Eq.  {B.  1). 

Two  particles  i  and  j  can  be  their  mutual  nearest 

neighbors,  if  there  is  no  other  particle  in  the  sp".ce  domain 

enclosed  by  two  spheres,  one  with  its  center  at  r.  and  with 

radius  r. .  and  the  other  with  its  center  at  r.  and  with  the  same 

J 

3 

radius.  See  Fig.  F.  1.  The  volume  of  the  domain  is  (9n/4)r^  . 
Kence 

W.j  =  exp[-  r..3  n]  (  B.  7) 

is  the  probability  that  particle  i  and  particle  j  are  their  mutual 
nearest  neighbors.  It  is  noted  that  the  above  consideration  is 
simply  a  matter  of  geometry. 
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APPENDIX  C 


Distribution  Functions  Decomposed  According  to  Inter -Particle  Distances 

The  distribution  function  F^(XjX2X3;t)  of  three  particles 
1,  2,  and  3  is  divided  into  parts  as  follows: 


f(3)  ^  F(3){X1X2X3;t) 

F(  1/2/3)  =  (1-W12)(1-W13(1-W23)  F(3) 

F(  12/3)  =  W1Z(1-W13)(1-W23)F<3) 

F(13/2)  =  (1-W12)W13(1-W23)F(3) 

F(l/23)  =  (1-W12)(1-W13)W23F(3) 

F(12/13)  =  W12W13(]-W23)F(3) 

F(13/2J)  =  (1“w12)W13W23F^^ 

F(  12/  13/23)  =  W12W15W23F(3) 

where,  for  example,  is  the  probability  o'f  particles  1  and  2 

being  their  mutual  nearest  neighbors.  The  sum  of  the  right  hand 
side  members  is  shown  to  be  F^  itself.  If  it  is  considered  that 
a  particle  has  only  one  nearest  neighbor  at  each  moment  of  time, 
we  may  have 


W12W23  ~  W23W12  =  W13W23  =  ° 


(C.  1) 


(C.2) 
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and  equations  tC .  i)  yield 


F(l/2/3)  =  (3  -  W12‘W13‘W23)f(3) 

F(i  z/y  =  w12f(3) 

F(13/2)  =  W13F(3)  (C,3) 

F(l/23)  =  W23F(3) 
and  it  holds  that 

F(l/2/3)  +  F(iZ/3)  +  F(13/Z)  +  F(l/23)  =  F(3) 

(C.4) 

Those  relations  given  by  \C.  3)  are  illustrated  in  Fig.  C.  1. 

If  we  consider  test  particle  1  and  four  field  particles  2,  3, 

4,  5,  the  possible  relations  among  them  are  illustrated  as  in  Fig.  C.2. 
By  taking  into  consideration  relations  similar  to  (C.2)  we  obtain 


F(5)  =  F(5)(X1X2X3X4X5;t) 
FO/2/3/4/5)  =  <1-W12)  (  1  -  Wu)  -  -  -  (1  •  W45)F 


(5) 


:U-  T.  W  +  E  WiiWk£. 

l<j  i  <j<k<4  J  £ 


(  1  -  Wl2  -  Wl3  -  W14  -  w15  -  W23  -  W24 


6.. 


*W25'W34'W35  ■  W45+  W12W34  +  W12  W35  +  W12W4S 


+  ws  w45)f 


(5) 


F(12/3/4/5)  =  W12(l-W34  -  W35  -  W45)  F 


(5) 


F(12/34/5)  =Wl2W34F 


(5) 


(C.5) 


etc. 


Finally  it  holds  that 

F(  1/2/3/4/5)  +  F(12/3/4/5) + - 

=  F^  (C.  6) 
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2.  T.  Koga,  Phys.  Fluids,  3,  454(1963);  Bull.  Am.  Phys. 

Soc.  11,  554(1966):  paper  presented  at  the  annual  meeting  of  the 

Division  of  Pla' ma  Physics,  APS,  1966. 

(N) 

3.  The  reason  for  giving  D  by  (2.  2)  has  been  discussed 

in  detail  in  the  second  paper  mentioned  in  footnote  2.  The  reader 

(N) 

might  th:  ->.k  that,  as  long  as  D'  '  is  given  by  (2.2),  there  is  no 

(N) 

need  of  consideration  of  Eq.  (2.  1)  since  D  given  by  (2.  2) 

is  a  solution  of  Eq.  (2.  1).  In  spite  of  its  simple  appearance, 

(N)  * 

D  given  by  (2.  2)  is  not  a  simple  function  of  t.  X  .  becomes 

known  only  by  solving  simultaneously  the  equations  of  motion  of 

the  N  particles;  these  equations  are  the  characteristic  equations 

of  Eq.  (2.  1),  and  the  sol.  nion  is  out  of  the  scope  of  our  mathematical 

ability. 


68 


4.  We  often  use  the  pres  nt  model  for  simulating  the  electron 
system  in  a  plasma.  The  simulation  is  feasible  when  our 
interest  is  in  the  interaction  between  a  plasma  and  an  external 
electro-magnetic  field  of  a  macroscopic  scale.  However, 

the  simulation  is  not  feasible  if  we  are  interested  in  effects 
of  inter-particle  interactions. 

5.  An  electron-electron  interaction,  for  example,  means  an 
interaction  of  two  electrons  which  are  their  mutual  nearest 
neighboring  particles. 

6.  See  footnote  1. 

7.  C.  Oberman,  A.  Ron,  and  J.  Dawson,  Phys.  Fluids,  5,  1514(196^), 
investigated  the  effect  of  an  external  force  of  a  short  time  scale, 

as  based  on  the  BBGKY  hierarchy,  but  the  assumption  of  time 
scales  necessary  for  the  hierarchy  was  not  properly  considered 
there. 

8.  Th;  problem  was  first  discussed  by  Hertz,  according 
to  S.  Chandrasekhar,  Rev.  Mod.  Phys.  15,  1  (1943). 


CAPTIONS  OF  FIGURES 


Fig.  3.  1.  A  field  particle  2  interacts  with  the  test  particle  1 
at  close  distances  with  time  scale  j\  after  the  interaction  particle 
2  interacts  with  part’~le  3,  particle  4,  and  so  forth.  The  time 
scale  of  those  close  interactions  is  the  same  as  T.  The 
space  domain  in  which  particle  2  is  localized  over  a  longer  period 
of  time  is  indicated  with  shadow.  The  force  exerted  on  particle 
1  by  particle  2  is  classified  as  follows:  1)  The  force  during  the 
close  interaction  particle  2  with  particle  1;  2)  the  force  during 

the  close  interactions  of  particle  2  with  particles  3,  4,  etc;  3) 
the  averaged  force  during  the  existence  of  particle  2  in  the 
shadowed  domain 

Fig.  4.  1.  In  A,  X ^  is  the  precise  trajectory  of  particle  i 

i 

under  the  influence  of  all  the  field  particles;  X.^  is  the  trajectory 
of  particle  i  when  the  influence  of  the  field  particles  is  ignored; 

II 

X.  is  the  trajectory  of  particle  i  interacting  with  particle  j  only. 

"...  *  (1) 

X.  is  much  closer  to  X^  .  In  B,  F^'  '  at  t  is  shown  by  a  6-function; 

<F.^>  is  the  result  of  averaging  from  t  to  t  +  t.  along  a 

i  / 1  \ 

set  of  trajectories  of  the  same  class  as  X.  ;  «F. '  '»  is  the 

n\ 

result  of  averaging  F^  ,  from  t  to  t  +  x,  along  a  set  of 

II 

trajectories  of  the  sam«  class  of  X^  . 
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1 


Fig.  6.  1.  Ion  I,  with  a  large  charge  and  a  large  mass,  interact 
with  a  similar  ion  2.  Many  electons,  depicted  with  small  dots, 
interfer  the  ion-ion  interaction. 

Fig.  B.  1  The  volume  of  the  space  domain  closed  with  two 

spheres  with  radius  r..  is  -2ll_  r.  ,3. 

K  ij  4  ij 

Fig.  C.  1.  Particle  1  is  the  test  particle;  particle  2  and  particle 
3  represent  the  field  particles.  There  are  four  typical  relations 
among  them:  a)  the  threes  are  mutually  remote;  b)  particle  i 
and  particle  2  are  their  mutual  nearest  neighbors;  c)  particle  1 
and  particle  3  are  their  mutual  nearest  neighbors;  d)  particle  2 
and  particle  3  are  their  mutual  nearest  neighbors. 

Fig.  C.  2.  Particle  1  is  the  test  particle;  particles  2,  3,4,  5, 
repiesent  the  field  particles ;among  them  particles  2  and  3  are 
ions  and  particles  4  and  5  are  electrons.  There  are  26  relations 
among  them;  full  and  short  lines  indicate  nearest  neighborhood 
and  long  and  dotted  lines  remote  relations. 
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FIG.  C.  1 
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